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ABSTRACT 

Data on low x physics from ZEUS & HI are presented and their interpretation 
discussed. The focus is on the increasing hardness of the energy dependence of 
inclusive 7*p scattering and certain diffractive processes as the transverse size of the 
probe decreases. 
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1 Introduction 



HERA remains a unique facility, the only e ± p collider. With 27.5 GeV beams on 
920 GeV protons the centre of mass energy is 318 GeV which also gives the upper 
bound on momentum transfer and thus the scale to which proton structure can be 
probed. Apart from extending the study of deep inelastic scattering and partonic 
physics, HERA also allows 7*p interactions to be study over a wide range of photon 
virtuality down to almost real photons. It provides a laboratory in which the high 
energy behaviour of cross-sections may be studied as a function of the transverse 
size of the projectile. This is the realm of 'low x' physics. This brief survey covers: 
data from HI and ZEUS on inclusive and diffractive processes: their analysis using 
perturbative QCD and other models, particularly colour diople models; the evidence 
or otherwise for universality and gluon saturation at low x. 



2 Formalism and phase space 

Inclusive electron-proton scattering at HERA ep — > eX is shown in Fig. [1] - LH. 
At a fixed centre of mass energy, -^/s, where s = (k + p) 2 ) the process is described 
by two kinematic variables, Q 2 , the four-momentum transfer squared and Bjorken 
x where 

Q 2 = -q 2 = ~{k - k') 2 : x = Q 2 /(2p.q). (1) 

The inelasticity, y = (p-q)/(p-k), the fractional energy loss of the electron in the 
proton rest frame, is related to x and Q 2 by Q 2 = sxy. The primary measured 
quantity is the double differential cross-section which, for Q 2 « M|, may be 
written in terms of two structure functions 

d 2 a 2ira 2 



dxdQ 2 xQ 4 



'l + (l-yf)F 2 (x,Q 2 )-y 2 F L (x,Q 2 ) . (2) 



The contribution from y 2 Fj J ^ the longitudinal structure function, is small and will 
be ignored in most of what follows. The inclusive scattering process may also be 
considered in terms of the total cross-section for 7*p scattering. At small x 

Air 2 n 

a%(W 2 ,Q 2 ) = —F 2 (x,Q 2 ), and W 2 « Q 2 /x, (3) 

where W 2 = (p + q) 2 is the centre of mass energy squared of 7*p system. So the 
structure function at small x gives cC? p at high CM energies. 

The x, Q 2 region in which inclusive measurements have been made is shown in 
Fig. |T| - RH. The strong correlation between x and Q 2 follows from the constraint 
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Figure 1: LH: Inclusive ep scattering and momenta. RH: Kinematic Plane for 
inclusive ep scattering up to HERA energies. 



Q 2 = sxy and the kinematic limit at HERA is given by the line y — 1. Fixed 
target measurements are in the lower right-hand quadrant of the plot. The region 
of interest for this talk is roughly x < 0.01, and Q 2 < 250 GeV 2 . 

3 Contexts 

The study of deep inelastic lepton scattering (DIS) and the development of QCD 
have been intimately related and this continues in the low x region opened up by 
HERA.0 The success of pQCD for hard processes relies on factorisation theorems 
which, in the case of DIS, enables the structure function F 2 to be written in the 
form 

F 2 (x,Q 2 ) = Y / e 2 f xq f (x,Q 2 ), with 9 f^ ] ~ g <U ® P «> ( 4 ) 

where qj is the momentum density function for quark of flavour / and charge e/, 
Pij are the QCD splitting functions and the sums run over q and q. At low x F 2 is 
largely given by the flavour singlet qq sea contribution and the evolution of this term 
is coupled to that of the gluon. The behaviour of both is dominated by the singular 
^^Much more detail and references to original papers may be found in review article of ref. [0. 
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behaviour of the gluon splitting function P gg ~ l/x, giving rise to an increase in the 
gluon density, the qq sea and hence F 2 at low x. 

The DGLAP equations (Eq. which have been calculated to next-to-leading order 
(NLO) in as and partially to next-to-next-to- leading (NNLO), represent a particu- 
lar, physically motivated, choice of how higher order corrections are summed. The 
most general behaviour of the splitting functions may be represented by double sum 
over terms of the form (In Q 2 ) n (ln l/x) m and there are other possible choices of sum- 
mation that may be relevant at low x. One of earliest results of pQCD as applied 
to DIS was derived by taking the leading log terms in both Q 2 and giving 

F 2 (x, Q 2 ) ~ exp ^(12a s /n)\n(l/x)\n(Q 2 /Q 2 ). (5) 

Another approach, taken by Balitsky et al (BFKL), involves summing the ln(l/x) 
terms and gives the striking prediction that F 2 ~ x~ x , where A ~ 0.5 at leading 
order and fixed as- For running as and NLO, the corrections to this result are 
large and are still subject to much debate. Despite the uncertainty surrounding the 
exact prediction, the BFKL approach has been of seminal importance in low x QCD 
dynamics. The DGLAP and BFKL approaches represent two 'extreme' choices in 
how the transverse momenta of the radiated gluons are ordered. All calculations 
agree that gluon dynamics dominates at low x and that the gluon density will 
increase as x decreases. At some point the gluon density will be so large that 
non-linear gluon recombination effects will need to be taken into account. A simple 
estimate suggests that such effects may occur for Q 2 values for which the gluon-gluon 
cross-section times the gluon density is of order the proton size: 

q 2 xg(x, Q ) ~ ttR , (6) 

where R is the proton radius and xg the gluon density, thus slowing the rise of F 2 . 
Hadronic total cross-sections involve soft physics which cannot be calculated per- 
turbatively. Regge theory provides a framework for describing the high energy be- 
haviour of total cross-sections. One expects crj^ drons ~ A ■ W 2 ^ ap ^ > ~ 1 \ where A is a 
process dependent constant and «p(0), the intercept of the 'Pomeron trajectory'^, is 
process independent. Although the value of «p(0) is not given, the model has been 
successfully applied to a wide range of hadron-hadron and real photoproduction 
data, giving a fitted value of 1.093(2) 0. If these ideas also apply to 7*p scattering, 
using Eq. [5], one might expect F 2 ~ a; 1 - ^ ) at small x. Using the hadronic value 

2 In principle ap is determined by the exchange of states with vacuum quantum numbers in the 
crossed channel. 
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of ap(0) predicts a rather gradual rise F 2 ~ x~ om , with the exponent independent 
of Q 2 . The BFKL calculation may be viewed as an attempt to calculate ctp(O) per- 
turbatively. Quasi-elastic scattering processes, such as diffraction, should also be 
dominated by the exchange of vacuum quantum numbers and the Pomeron, which 
is why they are important in unravelling the secrets of low x dynamics. 



4 F 2 at low x 

Fig. H shows F 2 data from HI, ZEUS and the NMC fixed target experiments, plotted 
as a function of x in bins of Q 2 . The RH plot, for Q 2 = 15 GeV 2 well into the deep 
inelastic range, shows F 2 rising steeply as x decreases. The LH plots show bins with 
Q 2 decreasing from 3.5 to 0.35 GeV 2 . Here the rise is still evident at the larger Q 2 
values, but diminishes rapidly as Q 2 decreases below ~ 1 GeV 2 . From Eq. |3| the 
rapid rise of F 2 at small x implies that a*°' p (W 2 , Q 2 ) increases more rapidly with 
W 2 as Q 2 increases, which is not expected in the hadronic Regge framework. 
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Figure 2: F 2 data from HI, ZEUS & NMC as a function of x in bins of Q 2 . 



HI have made a model independent study of the rise of F 2 ||. Using data with 
x < 0.01, the logarithmic slope with respect to lnx is calculated using 

<91nF 2 



HQ' 



d\nx 



or F 2 (x,Q 2 ) = c(Q 2 )x 



-A(Q ) 



(7) 



Q- 



The data from HERA are now sufficiently extensive and precise for A to be measured 
in bins of Q 2 over a range of x values. At fixed Q 2 , A is found to be independent of 
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x for x < 0.01. Fig. ^ shows A and c as functions of Q 2 where the errors have been 
calculated using the full systematic and statistical errors of the F2 data including 
correlations. At low values of Q 2 < 1 GeV 2 , the value of A is not far above the value 
expected from the hadronic Regge parameterisation, whereas at larger values of Q 2 , 
A rises almost linearly with InQ 2 to a value around 0.3 at Q 2 ~ 100 GeV 2 . The 
results for c(Q 2 ) are roughly independent of Q 2 with a value of about 0.18. The 
data show no indication that the rate of rise is starting to moderate at larger Q 2 , as 
might be expected from gluon saturation. Similar results for the logarithmic slope 
have been found in a preliminary analysis of ZEUS data 0]. 
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Figure 3: The logarithmic slope and constant versus Q 2 from the Fi{x,Q 2 



c(Q 2 )x 



-HQ 2 ) 



fit to data with x < 0.01. The F 2 data used is mostly from HI with 



some very low Q points from ZEUS. 



Although the results from the A analysis support the hypothesis of gluon dominance 
at low x, it is not possible from this study alone to say which of the pQCD approaches 
fits best. More information may be gained by applying the models discussed in Sec. |3| 
to the F 2 data over the full range of HERA and fixed target measurements. The 
most detailed fits are those using NLO DGLAP to determine the nucleon parton 
momentum densities. More details are given in the talk by Milstead ||, but the 
quality of the fits may be judged from Fig. which shows that NLO DGLAP can 
describe F2 for Q 2 values larger than 1.5 — 2 GeV 2 and all x values. For other 
pQCD approaches, fits to the low x F 2 data of comparable quality are obtained 
but the predictions for other observables, for example Fl, differ quite considerably. 
Although Fl has been extracted in the HERA kinematic region |J , the data do not 
yet have the precision to distinguish between the models. 
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5 Universality at low xl 



The F 2 data show gluon and qq sea dominance at low x. As the influence of the 
valence quarks is also small, the possibility arises that this behaviour could, in some 
sense, be 'universal'. To test this idea the structure function of another hadron must 
be extracted. Both ZEUS and HI have published data on the process ep — > enX 
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Figure 4: ep — > enX at HERA: (a) diagram for pion exchange; (b) the pion form 
factor extracted from the data at small t. From ffij. 



The neutron emerges at a very small angle with respect to the proton beam 
direction and is measured in a dedicated neutron counter at a distance of the order 
of 100 m downstream from the main detectors. Knowing the energy and angle of the 
neutron allows the momentum transfer squared at the pn vertex, t, to be calculated. 
At very small values of t, the process is dominated by charged pion exchange as 
shown in the LH diagram of Fig. |j. In this region of phase space the cross-section 
for 7*p — > nX may be written as a convolution of the flux of pions in the proton f n / p 
with the pion structure function F£ . The largest uncertainty in the measurement 
of F% is the pion flux. Current models differ by up to a factor of two, but this does 
not affect the x dependence which is shown in the RH plots of Fig. |] - in this case 
f n / p has been estimated using the additive quark model. The data show very clearly 
that F£ rises strongly for x values below 0.01. Thus at a qualitative level the data 
support the idea of a universal behaviour at small x. 



7 



6 Diffraction 



The diffractive process ep — *> epX is shown in Fig. |5|, 
together with the additional variables needed to de- 
scribe it beyond x, Q 2 and W 2 already defined: xp ~ 
1 —p' z /p z , the fractional longitudinal momentum loss 
of the proton; Mx the invariant mass of the diffractive 
final state at the 7* vertex; (3 = Q 2 /(2q.(p — p')) is 
the equivalent of Bjorken x in the fully inclusive case. 
Here the focus is on the W 2 and Q 2 dependence of two 
classes of diffractive processes: vector meson produc- 
tion for which Mx = My; inclusive diffraction with 



Figure 5: Diagram and ad- 
ditional kinematic variables 
for diffraction at HERA. 
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Mx > 3 GeV. Diffractive events are identified either 

by direct measurement of the scattered proton using a 'leading proton spectrometer' 
placed very close to the beam line at a large distance from the primary interaction 
or by the presence of a large rapidity gap in the main detector between the proton 
direction and the first energy deposits from the particles making up the Mx system. 
As already mentioned, diffraction is a quasi-elastic process and in many models the 
underlying physics is closely related to that of the elastic process. In particular in 
the Regge approach, the high energy behaviour of the total diffractive cross-sections 
should be controlled by Pomeron exchange. 

6.1 Vector Mesons 

First consider the data shown in the LH plot of Fig. |^, (Ttot{lP ~^ Vp) f° r rea l 
photoproduction of vector mesons and (TtotilP ~ * X) and for W > 10 GeV the 
data are fit to a power law dependence a ~ W . For the light vector mesons 
(p°,o>,0) 5 ~ 0.22, not inconsistent with the value expected from Regge theory. 
For 7£> — > J/ipp the energy dependence is much steeper, giving 5 ~ 0.8. The RH 
plot shows 5 from the energy dependence of j*p — > p°p in fixed Q 2 bins from to 
27 GeV 2 . Although the errors are large, there is a clear tendency for the energy 
dependence to steepen as Q 2 increases. In both cases the trend is clear, when there 
is a hard scale present, either a large My or a large Q 2 , then the cross-section rises 
more quickly with energy than expected from soft hadronic physics. 

6.2 Inclusive Diffraction 

Fig. [7| - LH shows a ZEUS measurement of the ratio of the diffractive cross-section 
to the total YP cross-section for two values of M x and 8 values of Q 2 between 0.27 
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Figure 6: LH plot, energy dependence of quasi- elastic vector meson photoproduction; 
RH plot, exponent of the energy dependence of cr(j*P ~ > P°p) vs Q 2 ■ 

and 60 GeV 2 as a function of W ||. Assuming Pomeron dominance, Regge theory 
predicts that the energy dependence of the ratio is given by (iy 2 ) 2ap " ap(0)_1 where 
ctp is the value of the Pomeron trajectory evaluated at the mean t of the diffractive 
data, for the hadronic Pomeron this gives an expectation of jy - 19 . Fitting the 
ratio data at fixed Mx and Q 2 to a form IV s gives the results 5 = 0.24 ± 0.07 for 
data with Q 2 < 1 GeV 2 and 5 = 0.00 ± 0.03 for data with Q 2 > 1 GeV 2 . Thus, 
once again, at Q 2 close to zero the data follows the expectation derived from soft 
hadronic physics, whereas in the deep inelastic region this is not the case and the 
diffractive cross-section (at fixed Mx and Q 2 ) follows the energy dependence of cr*°* . 



7 Colour Dipole Models 

Many models have been proposed to describe the transition from the soft hadronic- 
like energy dependence at Q 2 m to the harder behaviour seen at larger Q 2 . A 
particularly appealing approach is that offered by colour dipole models. In the 
rest frame of the proton the virtual photon splits into a qq pair a long time (or 
equivalently a large distance) before the interaction. The qq pair, characterised by a 
transverse size r and sharing the longitudinal momentum of the proton in the ratio 
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(1 — z) : z, then interacts with the proton with the cross-section o~ qq (x,r) giving 

a^ p {W\Q 2 ) = J d 2 rdz%,( r ,z,Q 2 )a qq (x,r)^ r (r,z,Q 2 ), (8) 
where \l/ 7 * is the known wave function for 7* — > gg. The dipole cross-section has 



to be modelled and a number of formulations have been proposed (see |L0| and 
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Figure 7: LH plots: Ratio of inclusive diffraction to a °* versus W at fixed M x and 



Q 2 . RH plot: The colour dipole model of Golec-Biernat et al.fit to F 2 data in bins 
of constant y. 



references therein). Here the version of Golec-Biernat & Wiisthoff (GBW) [O] is 
followed, partly for its simplicity, but also because it builds in the constraint of 
saturation. In their approach, the dipole cross-section is given by 



qq\ 



(x,r) = cr 1 - exp(-r 2 / ARl(x)) , R {x) = (x/x ) 



a/2 



(9) 



where a"o,Xo,A are three parameters to be determined from data. The interesting 
features of Eq. || are: for small dipole sizes r << i? , a qq oc r 2 x~ x ; for large dipoles 
r >> Rq, o qq oc cr and that the parameter R setting the scale depends on x. This 
last feature means that the approach to the 'saturation limit' occurs for smaller 
dipole sizes as x decreases. The parameters of the model are fixed by fitting the 
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HERA low x F 2 data. Another very compelling feature of the dipole approach is 
that diffraction may also be described. The diffractive cross-section is given by 



Diffr. 



dt 



167T 

t=0 



J d 2 rdz\y(r,z)\ 2 a 2 qq (x,r). (10) 



Thus no new parameters are required beyond those needed to describe the total 7*p 
cross-section. The successes of this simple model are a reasonable fit to the low x 
and low Q 2 (including the transition to the non pQCD region - see Fig. |7|-RH) and 
the energy independence of the ratio a^*p * r '/a!ft p [fl^| . An obvious deficiency of the 
model is that it does not contain any mechanism for pQCD Q 2 evolution of the 
structure function. An attempt to remedy this has been made by Bartels et al. |L3| . 
the dipole cross-section for small r is modified and related to the gluon density by 



7T 2 

<T qg (r,x) « —r 2 a s xg(x,fi 2 ), (11) 

where the scale fi 2 ~ C/r 2 and C is a parameter. The modified model now has 
5 parameters to be determined by the F 2 data. The fit to the higher Q 2 data 
is improved and \{Q 2 ) (Eq. 0) is well described for all Q 2 , whereas for the non- 
evolving model A falls below the slope data for Q 2 > 10 GeV 2 . The success of the 
GBW dipole model cannot be taken to imply that saturation is required by the 
HERA data as other models without saturation, e.g. the dipole model of Forshaw 
et al. |TIJ , give as good a representation of the data. 



8 Summary 

HERA has provided high precision data on the proton structure function F 2 over a 
wide range of x and Q 2 . In the low x region the bulk of the data is now systematics 
limit (typical uncertainties ~ 2%). HERA has also provided a wide range of mea- 
surements on hard diffractive scattering and quasi-elastic vector meson production. 
The striking feature of F 2 at low x - the strong rise - is mirrored in diffractive 
processes when there is a hard scale. The first measurement of the pion structure 
function at low x shows that it too rises strongly as x decreases, which hints at uni- 
versality in low x dynamics. The behaviour of low x F 2 , or equivalently high energy 
<r*°p, is dominated by gluon dynamics. Although some of the important calculations 
were completed well before HERA started operations, there is no doubt that HERA 
has opened up new avenues in strong interaction physics. Particularly the determi- 
nation of the gluon density at low x, the refinement of high density perturbative 
gluon dynamics and the deepening of the relationship between diffractive scattering 
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and the physics underlying the rise of total cross-sections with energy The HERA 
measurements and the related theoretical developments provide essential input for 
Run II at the Tevatron, RHICQ and the LHC. 
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